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ABSTRACT 

It is shown that in magnetohydrodynamics, the shapes of the magnetic 
energy and magnetic helicity spectra can influence whether the magnetic helicity 
resistive dissipation time scale is longer, equal to, or shorter than that of 
magnetic energy. The calculations highlight that magnetic helicity need not 
always dissipate significantly more slowly than magnetic energy. While this 
may be implicitly understood in detailed studies of magnetic helicity in plasma 
devices, it is not sufficiently emphasized in elementary discussions comparing 
magnetic helicity vs. magnetic energy evolution. Measuring magnetic energy 
and helicity decay times may provide a posteriori insight into the initial spectra. 

PACS codes: 52.30.Cv, 95.30. Qd, 52.65Kj, 52.55 -s, 96.60.Hv, 98.62Mw 



1. Introduction 



Magnetic helicity, (defined by / A • Bd 3 x, where A is the vector potential, B is 
the magnetic field and / d 3 x is a volume integral) is a measure of the linkage and 
twist of magnetic field lines, and is a conserved quantity for a closed system in ideal 
magnetohydrodynamics [1, 2]. Magnetic helicity is often more strongly conserved than 
magnetic energy in the presence of resistive dissipation [2-6]. The concept of a slower decay 
of magnetic helicity than magnetic energy forms the basis for Taylor relaxation theory [2-6], 
which describes the end state to which magnetically dominated magnetohydrodynamic 
configurations relax. The theory is well known in the context of laboratory pinches and 
has also been applied to the solar corona [7] and astrophysical jets [8]. Taylor theory 
predicts the relaxed state by minimizing the magnetic energy subject to the constraint that 
magnetic helicity is conserved. The result is a force-free configuration with the scale of the 
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field reaching the largest gradient scale available subject to boundary conditions. Some 
systems (e.g. Reverse Field Pinches (RFP)) transit toward and away from the Taylor state 
[e.g. 5], when driven. 

The relation between magnetic helicity and magnetic energy decay rates has been 
studied in [3] but without focusing on the role that energy and magnetic helicity spectra 
play in their relative magnitudes. Here I specifically show how the relative shapes of 
the magnetic and magnetic helicity spectra can influence the ratio of magnetic energy to 
magnetic helicity decay time scales. I derive the spectral shapes of magnetic helicity and 
magnetic energy for which the magnetic energy decays faster, equal, or more slowly than 
magnetic helicity. The purpose of these simple calculations is to illustrate very simply how 
the commonly intuited slow decay of magnetic helicity compared to magnetic energy is 
influenced by the spectral shapes of these quantities. 



2. Magnetic Helicity vs. Magnetic Energy Decay Time Scales 

In what follows, I will assume that surface terms are small compared to other terms 
involved in the time evolution of the the magnetic energy, magnetic helicity, and kinetic 
energies of the system. This would be the case if injection or loss of these quantities 
through the boundaries proceeds on time scales which are long compared to that of the local 
dissipation. It should be noted that this is not always the case [5, 9], but the assumption 
serves to simply illustrate the key concepts for the present purposes, and facilitates 
straightforward testing with spectral code numerical simulations. The time-evolution 
equations of magnetic helicity (e.g. [2]), magnetic energy and fluid kinetic energy for 
incompressible, resistive, viscous MHD can then be written 

d t H = d t (A-B) = -2X(3-B), (1) 

d t M = d t (B 2 ) = -2A((VB) 2 ) - 2(v • (J x B)>, (2) 

and 

d t K ee d t (v 2 ) = -2z/((Vv) 2 ) + 2(v • (J x B)>, (3) 

where H is the magnetic helicity, M is twice the magnetic energy, K is twice the kinetic 
energy, A is the magnetic diffusivity, v is the viscosity, v is the velocity, and B is the 
magnetic field in Alfven units. The vector potential A satisfies B = VxA, and the current 
density is normalized such that J = VxB. The brackets indicate a volume average. 

In general, the velocity should be included dynamically, as is done in [10]. Here I 
consider a case in which the velocity term can be approximately ignored in the energy 
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dynamics. Deriving a condition for this requires an analysis of the spectra, as discussed in 
the appendix. Note that even if such a condition were not satisfied, the magnetic energy loss 
rates and subsequent comparisons to those of magnetic helicity would still be of interest: 
they would specifically represent the resistive dissipation of magnetic energy, which would 
not necessarily be approximately equal to the total magnetic energy change rate. 

Now write M = M k dk with 

M k = J \B(k)\ 2 k 2 dQ k = J k 2 \A(k)\ 2 k 2 dQ k oc k~ q , (4) 

where the tilde indicates the Fourier transform, B(k) = i\t x A(k), k is the wavenumber 
magnitude, and Q k is the solid angle in wavenumber space. The second relation in (4) 
follows in the Coulomb gauge V ■ A = 0, and the assumed power-law form of M k takes q to 
be independent of k between ki and k\ » ki, the smallest and largest wave numbers of 
the system (e.g. the inverse of the system size, and the resistive wavenumber respectively). 
Similarly, we write the magnetic helicity spectrum 

H k = l -J [A(k) • B*(k) + A*(k) • B(k)]A; 2 ^ fc = M k f(k)/k, (5) 

where 

f(k) oc k~ s (6) 

is the assumed form of the fractional magnetic helicity as a function of k, with s independent 
of k, and * indicates complex conjugate. In the Coulomb gauge, / is the fractional helicity 
where |/| < 1 (realizability condition [11]), and |/| = 1 is the maximal helicity case. (Ref. 
[12] provides further discussion in the context of numerical simulations.) We also write the 
current helicity spectrum C k = kM k f(k). The above definitions imply 

(A • B) = H = f kx H k dk = [ kx f(k)M k k~ l dk, (7) 

Jk L Jk L 

(J • B) = C = / k 2 H k dk = / f(k)kM k dk, (8) 

J k L Jk L 

(B 2 ) — M — M k dk, (9) 

Jk L 

and 

((VB) 2 ) ~ / " k 2 M k dk. (10) 

Jk L 

These formulae can be used to calculate the ratio of the time scale of magnetic helicity 
dissipation to that of magnetic energy dissipation. In general, f(k) could have both positive 
and negative signs, but for simplicity, we consider an initial condition in which f(k) is 



-4- 



positive for all k within the integration bounds. Justification comes from 3-D numerical 
simulations of closed turbulent systems which show that unless kinetic helicity is injected, 
the magnetic helicity basically maintains the same sign over the spectral range, though 
transfer between scales does occur [13, 14]. 

For the magnetic helicity dissipation time scale, we then have 

tf x f(k)M k k~ l dk , N 

TH = —Hd t H = H/2XC = JkL \ K) , 11 

2\tj{k)kM k dk ) 

while for the magnetic energy decay time scale (assuming resistive decay) we have 

r k \ jy[ ^ 

r u K -MI8 t M K -MH3.MU = 2 >^J Mtdk . (12) 

In Fig. 1, the ratio R = t h /t m is plotted as a function of q and s using (11) and (12), for 
four different values of kj^/kx- The plots show a significant parameter regime where R < 1. 



3. Discussion 

3.1. Understanding specific regimes analytically 

To better understand Fig. 1, it is useful to further study the ratio R analytically. 
Using (4) and (6) and assuming q ^ 1,3 and q + s ^ 0, 2, we have 

T H = ~ " —, 13 

2(q + s)X[k^tl 

and 

™ ~ 2A(l- ? )[*Mif <14) 

where the brackets denote taking the difference between the quantity inside evaluated at k\ 
and I now explore several different value ranges of q and s and compute t m and r H and 
combine the results to determine R = t h /t m . 

First consider th- For < q + s < 2, Eqn. (16) gives 

~ 2( 9 + s) A*f-'-»' <15) 
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while for q + s > 2, 

(2 - q - s 
2(q + s)Xk 2 L 

Now consider tm- For 1 < q < 3, Eqn. (14) gives 



2\(1-q)h 

while for g > 3, 



_ (3 ~ q)k l L ~ q h?) 



(3-g) 
2A(l-?)*£ 

and for < q < 1 



^ = :^ 2 , (is) 



(3-g) nq v 

™ = 2A(i-«z)*r (19) 

Combining (15) and (17) for < g + s < 2 and 1 < q < 3, we have 

R = Q(k x /k L ) 1+s , (20) 

where Q = | ^"^"^ [ . Combining (15) and (18) for < q + s < 2 and q > 3, we have 

R = Q{k\/ ki) q+s ~ 2 . (21) 

Combining (15) and (19) for < q + s < 2 and < q < 1, we have 

R = Q(k x /k L y +s . (22) 

Combining (16) and (17) for q + s > 2 and 1 < q < 3, we have 

i* = Q(k x /k L f-i. (23) 

Combining (16) and (18) for q + s > 2 and g > 3, we have 

i? = Q. (24) 

Finally, combining (16) and (19), we have for q + s > 2 and < q < 1 

i? = Q(h/k L ) 2 . (25) 

Three qualitatively different regimes are revealed by (20-25): I. R < 1: Eqns. (20) and 
(21) show cases in which the magnetic helicity can actually decay more quickly than the 
magnetic energy. Eqn. (20) shows that for s < —1 and < s + q < 2 and 1 < q < 3, 
R < 1 can arise. In addition, Eqn. (21) shows that for < q + s < 2 and q > 3, i? < 1 is 
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guaranteed. II. R ~ 1: Eqn. (24) shows that R is independent of the ratio of wave numbers 
when q + s > 2 and q > 3. For this regime, Q > 1 so that the magnetic helicity decays 
more slowly, but Q can be close to 1, so R ~ 1. III. R > 1: For the remaining cases of Eqn. 
(20) for s > — 1, and for all regimes of q, s when Eqns. (22), (23) and (25) are applicable, 
t~h > tm by a factor that depends on k\/ki- These three regimes can be seen within the 
continuous surfaces of Fig. 1, where R was calculated directly from (11) and (12). 

In regime III, magnetic energy decays faster than magnetic helicity, and thus will 
indeed relax toward the Taylor state [3-5]. As just derived, this regime arises when both the 
magnetic energy and fractional magnetic helicity decrease with increasing k, and only when 
the magnetic energy spectrum is not too steep; regime II rather than regime III applies for 
a wide range of s > 0, when q > 3. Note that in contrast to regimes II and III, regime I can 
only result from an initial spectrum in which the fractional magnetic helicity increases with 
decreasing wave number. The realizability condition (\f(k) < 1|) [11] then further implies 
the system is not maximally helical at all scales when case I applies. This in turn implies 
that there would be scales at which the kinetic energy can equal or dominate the magnetic 
energy. An example that could qualify for regime I would be one in which the kinetic energy 
and magnetic energy were equal at large scales, but the velocity falls off more rapidly than 
the magnetic energy at smaller scales. Whether or not a particular physical phenomenon 
produces such given spectrum is not addressed here. The main point is just to highlight 
that the slow resistive decay is of magnetic helicity is not guaranteed in all circumstances. 

A subtlety occurs if f(k) were to change rapidly in time (e.g. if its spectral index s 
changes with time). The estimated th could then be time dependent in such a way that the 
initial reduction of H by a factor of 1/e could occur in a different parameter regime from 
the subsequent time evolution of H. For example, H could decay to a value of 1/e faster 
than M reduces by the same fraction so that R < 1 initially, but the subsequent evolution 
to H = might be slower than M if the change in f(k) induces a regime change such that 
R > 1 at late times. 

In principle, the results above could be used to help diagnose the initial fractional 
magnetic helicity and magnetic energy spectra for a system in which the decay time scales 
of these quantities can be measured. The ratio th/tm can then be compared with unity 
and the regime I, II or III revealed. The results can certainly be tested with a 3-D numerical 
experiment in which magnetic helicity and energy spectra are imposed initially to be in 
either regime I, II, or III and the time scale for the decay of magnetic energy or magnetic 
helicity monitored, and a comparison of the resistive decay terms measured. 
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3.2. Relation to previous work 

It is useful to discuss the implications herein in the context of Ref. [3]. The results 
of the two approaches are compatible, but complementary. Using the Schwartz inequality, 
Ref. [3] derived an upper bound to the total helicity decay rate as a function of the 
magnetic energy and magnetic energy decay rate, which was then used to determine how 
much magnetic helicity AH would decay in the time scale that magnetic energy decays, 
or, for a given amount of helicity, how long it would take to decay this amount of helicity. 
No explicit dependence of these decay rates on the spectral indices was considered. In the 
specific example of a coronal loop considered in Ref. [3], the helicity decayed more slowly 
than magnetic energy, which was consistent with that example reflecting a more nearly 
maximal helicity case. Indeed as discussed above, th > tm in regime III. But above it was 
revealed that r H = t m in regime II and r H < r M in regime I. These latter two regimes were 
not explicitly considered in Ref. [3]. 

The existence of the different regimes suggests that the commonly stated concept that 
"magnetic helicity decays on a longer time scale than magnetic energy" is not universally 
applicable. However, the modified statement "magnetic helicity resistively decays on the 
same or on a longer time scale than magnetic energy" is fully justified when based on 
a comparison of H max /d t H (or H max / d t H max ) to M/d t M, where H max is the maximum 
possible magnetic helicity (i.e. corresponding to f(k) = 1 for all k) for a given total 
magnetic energy. In this case, region I does not apply, and only region II and III apply. 
Note however, that a given system is not guaranteed to have maximal helicity. 

A comparison of H max /d t H to M/d t M also follows from equation (15) of Ref. [3], the 
key equation relating the two decay rates using the Schwartz inequality. Using the present 
notation, the content of that equation is 

H max /d t H > (M/dtM) 1 / 2 , (26) 

which is fully compatible with, but complementary to, the spectral considerations above. It 
does not by itself directly reveal the role of the spectral shapes. 



4. Conclusions 

Assuming that boundary terms are small, the influence of the spectral shapes of 
magnetic helicity and magnetic energy on the resistive decay time scales of total magnetic 
helicity and magnetic energy was determined. Regimes exist for which the magnetic helicity 
decays more slowly, more quickly, or at the same rate as magnetic energy. The conventional 
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regime for which r H /r M > 1 applies to a wide range of spectra, but the existence of regimes 
for which t h /t m = 1 and t h /t m < 1, while implicitly understood in some treatments of 
magnetic helicity evolution in laboratory devices [e.g. 5], are often not emphasized when 
the basics of magnetic helicity vs. magnetic energy decay are introduced [6]. When the 
magnetic helicity is nearly maximal, th/tm > 1, with the equality applying only when the 
magnetic energy spectrum is steep. 

Because I have assumed the velocity term in (2) is small (see also the appendix) care 
must be taken when evaluating the implications for general systems incurring dynamical 
Taylor relaxation in which velocity fluctuations must be solved for dynamically. Indeed 
Eq. (2) reveals that the velocity term provides an additional way for magnetic energy to 
drain into kinetic energy and subsequently dissipate through viscosity. This can help justify 
the conventional intuition that magnetic helicity decays more slowly than magnetic energy 
under certain circumstances and it also plays an important dynamical role in driving a 
system toward and away from the Taylor state [5, 10]. However, since the velocity term can 
be of either sign, it may also induce a slower decay of magnetic energy depending on the 
situation. In the presence of dynamically significant velocity fluctuations, the comparison 
made of time scales in this paper can then be simply thought of as tjj/ (TM) res , where the 
denominator represents the magnetic energy loss time through resistive dissipation only. 



It is helpful to compare the relative magnitudes of the terms on the right of (2) to 
motivate ignoring the velocity term. To do so, we write the kinetic energy as 



where E k oc k~ qv . We first work on the second term on the right of (2). Noting that 
(v ■ (J x B)) < ((v 2 ) 1 / 2 /c 2 B 2 ). and using the spectral scalings for kinetic energy combined 
along with M k oc k~ q from (4), we obtain 



Appendix: On ignoring the velocity term in Eq. (2) 




(27) 



(JxB))<k L B 2 L v L (k L /k) 



2q+qy-5 
2 



(28) 



so that for q < 



5-qy 
2 



(J x B)> ~ k L B 2 L v L {k u /k L ) 



5-2q-qy 
2 



(29) 



while for q > 



5-<?y 
2 



(v • (J x B)) ~ k L B\v L . 



(30) 
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In the above equations, the viscous cutoff wavenumber k v has entered because I have 
assumed that the magnetic Prandtl number Ptm = v/\» f: In magnetically dominated 
systems such as coronae or laboratory devices, Ptm ~ (T/1.4 x 10 5 K)(n/10 16 cm~ 3 ) -1 >> I 
[15]. Thus the maximum wavenumber for which (v • (J x B)) remains finite is k v > k\, 
explaining the appearance of k v . 

For the first term on the right of (2), we again use oc k~ q as in (4). We also assume 
that the decay time at k\, the dissipation scale, is approximately equal to the dynamical 
time at that scale, equivalently, 

\hl~h x B kx , (31) 
remembering that B is in Alfven units. It follows from some straightforward algebra that 

A((VB) 2 ) ~ k L B\R^= k L Bl(k x /k L )^, (32) 

where the magnetic Reynolds number R M = B L jki)\. In this approach, the velocity term in 
(2) can be dropped when Eqn. (32) >> Eqn. (29) or Eqn. (30), depending on the regime. 
More specifically, this condition becomes 

— (k L /K) q ~^ (h/K)^ » 1 (33) 
vl 

when q < and 

^(WM^»1 (34) 

when q > ^y 11 . Note that the square root of the ratio of magnetic to kinetic energies which 
enters the above relations is that at the outer scale, not the total, so the total magnetic 
energy can be larger than the total kinetic energy even the two energies do not satisfy this 
ordering at any given scale. 

To see what the conditions above would imply for a specific example, consider Bl ~ vl 
and let qy = 5/3. Then = 5/3, so for q < 5/3, the condition (33) would always be 
satisfied in this case. It is common to find steeper magnetic energy spectra when compared 
to kinetic energy spectra in MHD turbulence [14, 16, 17, 18]. 
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Figure 1: R = th /tm vs. fractional helicty spectral index s, and energy spectral 
index q. The plots were obtained by taking the ratio of (11) to (12) for 4 different 
values of k\/k L . Clockwise from top left: h\jh L = 10, 10 2 , 10 4 , 10 3 . For s > and 
3 > q > 0, R » 1. This is the parameter regime III in the text, and corresponds 
to the commonly studied regime where magnetic helicity decays more slowly than 
magnetic energy. However, note that for small q and s < 0, R < 1, and for q > 3, 
R < 1. These correspond to regimes I and II discussed below (25). The flat-top 
surfaces on each of the plots are just an artifact of truncating the vertical axis below 
the peak values of R » 1 attained by the actual solution surfaces. 



